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Abstract 

Cumulative resource constraints can model scarce resources in scheduling prob- 
lems or a dimension in packing and cutting problems. In order to efficiently solve 
such problems with a constraint programming solver, it is important to have strong 
and fast propagators for cumulative resource constraints. One such propagator is the 
recently developed time-table-edge-finding propagator, which considers the current 
resource profile during the edge-finding propagation. Recently, lazy clause gener- 
ation solvers, i.e., constraint programming solvers incorporating nogood learning, 
have proved to be excellent at solving scheduling and cutting problems. For such 
solvers, concise and accurate explanations of the reasons for propagation are essential 
for strong nogood learning. In this paper, we develop the first explaining version 
of time-table-edge-finding propagation and show preliminary results on resource- 
constrained project scheduling problems from various standard benchmark suites. 
On the standard benchmark suite PSPLib, we were able to close one open instance 
and to improve the lower bound of about 60% of the remaining open instances. 
Moreover, 6 of those instances were closed. 



1. Introduction 

A cumulative resource constraint models the relationship between a scarce resource and 
activities requiring some part of the resource capacity for their execution. Resources can 
be workers, processors, water, electricity, or, even, a dimension in a packing and cutting 
problem. Due to its relevance in many industrial scheduling and placement problems, it 
is important to have strong and fast propagation techniques in constraint programming 
(Cp) solvers that detect inconsistencies early and remove many invalid values from the 
domains of the variables involved. Moreover, when using Cp solvers that incorporate 
"fine-grained" nogood learning it is also important that each inconsistency and each 
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Figure 1: Five activities with prece- 
dence relations. 
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Figure 2: A possible schedule of the 
activities. 



value removal from a domain is explained in such a way that the full strength of nogood 
learning is exploited. 

In this paper, we consider renewable resources, i.e., resources with a constant re- 
source capacity over time, and non-preemptive activities, i.e., whose execution cannot 
be interrupted, with fixed processing times and resource usages. In this work, we de- 
velop explanations for the time-table-edge-finding (TtEf) propagator [34] for use in lazy 
clause generation (Lcg) solvers [22| [9] . 

Example 1.1. Consider a simple cumulative resource scheduling problem. There are 5 
activities A, B, C, D, and E to be executed before time period 10. The activities have 
processing times 3, 3, 2, 4, and 1, respectively, with each activity requiring 2, 2, 3, 2, 
and 1 units of resource, respectively. There is a resource capacity of 4. Assume further 
that there are precedence constraints: activity B must finish before activity D begins, 
written B <C D, and similarly C <C E. Figure [T] shows the five activities and precedence 
relations, while Fig. [2] shows a possible schedule, where the start times are: 0, 0, 3, 5, 
and 5 respectively. 

In Cp solvers, a cumulative resource constraint can be modelled by a decomposition 
or, more successfully, by the global constraint cumulative [2]. Since the introduction 
of this global constraint, a great deal of research has investigated stronger and faster 
propagation techniques. These include time-table [2j, (extended) edge-finding [2"T| [55]. 
not-first/not-last [21, 25j, and energetic-reasoning propagation [HE]. Time-table propa- 
gation is usually superior for highly disjunctive problems, i.e., in which only some activi- 
ties can run concurrently, while (extended) edge- finding, not-first /not-last, and energetic 
reasoning are more appropriate for highly cumulative problems, i.e., in which many ac- 
tivities can run concurrently. [4] The reader is referred to [6] for a detailed comparison of 
these techniques. 

Vilim [3l] recently developed TtEf propagation which combines the time-table and 
(extended) edge-finding propagation in order to perform stronger propagation while hav- 
ing a low runtime overhead. Vilim [33] shows that on a range of highly disjunctive open 
resource-constrained project scheduling problems from the well-established benchmark 
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library PSPLib|J TtEf propagation can generate lower bounds on the project dead- 
line (makespan) that are superior to those found by previous methods. He uses a Cp 
solver without nogood learning. This result, and the success of Lcg on such prob- 
lems, motivated us to study whether an explaining version of this propagation yields an 
improvement in performance for LCG solvers. 

In general, nogood learning is a resolution step that infers redundant constraints, 
called nogoods, given an inconsistent solution state. These nogoods are permanently or 
temporarily added to the initial constraint system in order to reduce the search space 
and/or to guide the search. Moreover, they can be used to short circuit propagation. 
How this resolution step is performed is dependent on the underlying system. 

Lcg solvers employ a "fine-grained" nogood learning system that mimics the learning 
of modern Boolean satisfiability (Sat) solvers (see e.g. [SO])- In order to create a strong 
nogood, it is necessary that each inconsistency and value removal is explained concisely 
and in the most general way possible. For Lcg solvers, we have previously developed 
explanations for time-table and (extended) edge-finding propagation |27j . Moreover, for 
time-table propagation we have also considered the case when processing times, resource 
usages, and resource capacity are variable |24j . Explanations for the time-table propaga- 
tor were successfully applied on resource-constraint project scheduling problems [27| [29] 
and carpet cutting [28] where in both cases the state-of-the-art of exact solution methods 
were substantially improved. The explanations defined here are similar to the step-wise 
ones for the (extended) edge-finding propagation in [27], but there we do not consider 
the resource profile and are more complex. Moreover, the proposed explanations for 
edge-finding propagation in [27] has never been implemented. 

Explanations for the propagation of the cumulative constraint have also been pro- 
posed for the PaLM [3 [13] and SCIP [TJ H2] frameworks. In the PaLM framework, 
explanations are only considered for time-table propagation, while the SCIP framework 
additionally provides explanations for energetic reasoning propagation and a restricted 
version of edge-finding propagation. Neither framework consider bounds widening in or- 
der to generalise these explanations as we do in this paper. Other related works include 
(32], which presents explanations for different propagation techniques for problems only 
involving disjunctive resources, i.e., cumulative resources with unary resource capacity, 
and generalised nogoods [15] . A detailed comparison of explanations for the propagation 
of cumulative resource constraints in Lcg solvers can be found in |24j . 

In this paper we develop explanations for the TtEf cumulative propagator in Lcg 
solvers. The explaining TtEf propagation is then compared with the explaining time- 
table propagation from |27] in the LCG solver on RCPSP using the reengineered LCG 
solver [9J which was also used for the experiments presented in |27j . 



2. Cumulative Resource Scheduling 

In cumulative resource scheduling, a set of (non-preemptive) activities V and one cu- 
mulative resource with a (constant) resource capacity R is given where an activity i is 
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specified by its start time Si, its processing time pi, its resource usage rj, and its en- 
ergy ei := pi ■ Ti. In this paper we assume each Si is an integer variable and all others 
are assumed to be integer constants. Further, we define esti (ecti) and Isti (Icti) as the 
earliest and latest start (completion) time of i. 

In this setting, the cumulative resource scheduling problem is defined as a constraint 
satisfaction problem that is characterised by the set of activities V and a cumulative 
resource with resource capacity R. The goal is to find a solution that assigns values 
from the domain to the start time variables Si (i G V), so that the following conditions 
are satisfied. 

esti < Si < Isti, Vi € V 

n < R Vr 

lev-.reiSi^Si+Pi) 

where r ranges over the time periods considered. Note that this problem is NP-hard [5]. 

We shall tackle problems including cumulative resource scheduling using Cp with 
nogood learning. In a Cp solver, each variable Si, i 6 V has an initial domain of possible 
values D°(Si) which is initially [esti, Isti]. The solver maintains a current domain D 
for all variables. Cp search interleaves propagation with search. The constraints are 
represented by propagators that, given the current domain D, creates a new smaller 
domain D' by eliminating infeasible values. The current lower and upper bound of the 
domain D(Si) are denoted by lb(Si) and ub(Si), respectively. For more details on Cp 
see e.g. |23j. 

For a learning solver we also represent the domain of each variable Si using Boolean 
variables {S{ < v\,esti < v < Isti. These are used to track the reasons for propagation 
and generate nogoods. For more details see [22]. We use the notation jv < Si}, esti < 
v < Isti as shorthand for ->{Si < v—ij, and treat \v < SiJ, v < esti and \S% < v J , v > Isti 
as synonyms for true. Propagators in a learning solver must explain each reduction in 
the domain by building a clausal explanation using these Boolean variables. 

Optimisation problems are typically solved in Cp via branch and bound. Given an 
objective obj which is to be minimised, when a solution is found with objective value o, 
a new constraint obj < o is posted to enforce that we only look for better solutions in 
the subsequent search. 



3. TTEF Propagation 

In this section we develop explanations for TtEf propagation. For a more detailed 
description about TtEf propagation the reader is referred to |34j . 

TtEf propagation splits the treatment of activities into a fixed and free part. The 
former results from the activities' compulsory part whereas the latter is the remainder. 
The fixed part of an activity i is characterised by the length of its compulsory part pf T := 
max(0, ecti — Isti) and its fixed energy ef T := ■ pJ T ■ The free part has a processing 
time pf F := pi — pf T and a free energy of ef F := — ef F . Let V EF be the set of 
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Figure 3: A diagram illustrating an activity i when started at esti or Isti, and its possible 
range of start times, as well as the compulsory part CPi, and the fixed and free 
parts of the processing time. 



activities with a non-empty free part {i G V | pf > 0}. An illustration of this is shown 
in Figure [3j 

TtEf propagation reasons about the energy available from the resource and energy 
required for the execution of activities in specific time windows. The start and end times 
of these windows are determined by the earliest start and the latest completion times of 
activities i 6 V . These time windows [begin, end) are characterised by the so-called 
task intervals V EF (a,b) := {i 6 V EF \ est a < esti A Icti < lctb} where a, b G V EF , 
begin := est a , and end := lctb- 

It is not only the free energy of activities in the task interval V EF {a,b) that is 
considered, but also the energy resulting from the compulsory parts in the time win- 
dow [estb,lctb). This energy is defined by ttEn(a,b) := ttAfter[est a ] — ttAfter[lctb] 
where ttAfter[r] := £ t > T T,ieV:lsU<t<ec tl r i- 

Furthermore, we also consider activities i € V \ V EF (a, b) in which a portion of their 
free part must be run within the time window as described in |34j . Let lstf F be the 
latest start time of the free part of an activity, i.e., Icti — p EF ■ Then activity i's free 
part consumes at least • (Ictb — lstf F ) energy units in [est a ,lctb) if est a < esti and 
lst EF < lct b . We define the energy contributed by such activities by rsEn{a, b) := 
'Eiev\v^(a,b):est a <esu max(0, lct b - lst EF ). 

In summary, TtEf propagation considers three ways in which an activity i can 
contribute to energy consumption within a time window determined by a task inter- 
val V EF (a, b). First, the free parts that must fully be executed in the time window; 
second, the compulsory parts that must lies in the time window; and third, some free 
parts that must partially be run in the time window. Thus, the considered length of an 
activity i is 

Pl ieV EF (a,b) 

Pi(a, b) := I max(0, lctb - IsU) i V EF (a, b) A est a < esti 

k max(0, mm(lctb, ecti) — max(es£ a , Isti)) others 

The considered energy consumption is e«(a, b) := ri ■ Pi{a, b) in the time window. 
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3.1. Explanation for the TTEF Consistency Check 

The consistency check is one part of TtEf propagation that checks whether there is a 
resource overload in any task interval. 

Proposition 3.1 (Consistency Check). The cumulative resource scheduling problem is 
inconsistent if 

R ■ (lctb — est a ) — energy(a, b) < — > _L (1) 
where energy(a, b) := X^ev BF (a b) e f F + ttEn(a, b) + rsEn(a, b). 

This check can be done in 0(l 2 + n) runtime, where I = (V^l, if the resource profile 
is given. The corresponding algorithm is shown in Alg. [T] in App. |A} 

A naive explanation for a resource overload in the time window [est a ,lctb) only con- 
siders the current bounds on activities' start times S{. 

f\ lest, < Si} A [St < Istij -> J_ 

i£V:pi(a,b)>0 

However, we can easily generalise this explanation by only ensuring that at least Pi(a, b) 
time units are executed in the time window. This results in the following explanation. 

/\ \est a + Pi (a, b)-pi< S t j A [Si < lct b - Pi(a, b)j -> J_ 

ieV:p,(a,6)>0 

Note that this explanation expresses a resource overload with respect to energetic rea- 
soning propagation which is more general than TtEf. 

Let A := energy(a,b)—R-(lctb—est a ) — l. If A > then the resource overload has extra 
energy. We can use this extra energy to further generalise the explanation, by reducing 
the energy required to appear in the time window by up to A. For example, if fj < A 
then the lower and upper bound on Si can simultaneously be decreased and increased by 
a total amount in {1, 2, minQA/VjJ ,Pi(a, b))} units without resolving the overload. If 
rj • pi(a,b) < A then we can remove activity i completely from the explanation. In a 
greedy manner, we try to maximally widen the bounds of activities i where Pi(a, b) > 0, 
first considering activities with non-empty free parts. If Aj denotes the time units of the 
widening then it holds Pi(a, b) > Aj > and J2ieV-pi(a b)>o Aj • r» < A and we create the 
following explanation. 

/\ lest a + Pi (a, b)- Pi -Ai< St} A [Si < lct b - Pi(a, b) + Aj] ->• !_ 

i£V:pi(a,b)-Ai>0 

The last generalisation mechanism can be performed in different ways, e.g. we could 
widen the bounds of activities that were involved in many recent conflicts. Further study 
is required to identify which are the most appropriate. 



6 



3.2. Explanation for the TTEF Start Times Propagation 



Propagation on the lower and upper bounds of the start time variables Si are symmet- 
ric; Consequently we only present the case for the lower bounds' propagation. To prune 
the lower bound of an activity u, TtEf bounds propagation tentatively starts the ac- 
tivity u at its earliest start time est u and then checks whether that causes a resource 
overload in any time window [est a ,lctb) ({a,b} C V ). Thus, bounds propagation and 
its explanation are very similar to that of the consistency check. 

The work of [34J considers four positions of u relative to the time window: right (est a < 
est u < Ictf) < ect u ), inside (est a < est u < ect u < Ictb), through (est u < est a A Ictb < 
ect u ), and left (est u < est a < ect u < Ictb)- The first two of these positions correspond 
to edge-finding propagation and the last two to extended edge-finding propagation. We 
first consider only the right and inside positions, i.e., est a < est u . Note that a could be 
u. Then, 

rest(a, b, u) 



R ■ (Ictb — est a ) — energy(a, b, u) < 



<S U 



where energy(a, b, u) := energy(a, b) — e u (a, b) + r u • (min(Zct&, ect u ) — est u ) and 



rest(a, b, u) := energy(a, b, u) — (R — r u ) ■ (Ictb — est a ) 

- r u ■ (mm(lctb, ect u ) - est u ) . 

The first two terms in the sum of energy(a,b,u) gives the energy consumption of all 
considered activities except u, whereas the last term is the required energy of u if it is 
scheduled at est u in the time window [est a , Ictb)- The propagation, including explanation 
generation, can be performed in 0(l 2 + k-n) runtime, where I = \V \ and k the number 
of bounds' updates, if the resource profile is given. Moreover, TtEf propation does not 
necessarily consider each u G V , but those only that maximise min(e^ F ,r u ■ (Ictb ~ 
est a )) — r u ■ max(0, Ictb — ^ s tu F ) an d satisfy est a < est u . The corresponding algorithm is 
shown in Alg. [2] in App. [Aj 

A naive explanation for a lower bound update from est u to newLB := \rest(a, b, u)/r u ] 
with respect to the time window [est a ,lctb) additionally includes the previous and new 
lower bound on the left and right hand side of the implication, respectively, in comparison 
to the naive explanation for a resource overload. 

\est u < S U J A /\ [estt < Si} A [Si < IsUj \newLB < S u j 

ieV\{u}: Pi (a,b)>0 

As we discussed in the case of resource overload, we perform a similar generalisation 
for the activities in V \ {«}, and for u we decrease the lower bound on the left hand 
side as much as possible so that the same propagation holds when u is executed at that 
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Table 1: Specifications of the benchmark suites. 



suite 


sub-suites 


T^inst 


^fcact 


Pi 


#res 


notes 


AT [3] 


st27/st51/st103 


48 each 


25/49/101 


1-12 


6 each 




PSPLib Qi] 


J30 [T7]/j60/j90 


480 each 


30/60/90 


1-10 


4 each 






j120 


600 


30 


1-10 


4 




BL g] 


bl20/bl25 


20 each 


20/25 


1-6 


3 each 




Pack [8] 


55 


15-33 


1-19 


2-5 




KSD15_d [IH] 


480 


15 


1-250 


4 


based on j30 


Pack_d [18 




55 


15-33 


1-1138 


2-5 


based on 
Pack 



decreased lower bound. 

\est a + Ictb — new LB + 1 — p u < S U }A 

f\ {est a + pj (a, b)-pi< Si} A {Si < lct b - pi(a, b)j 

ieV\{u}: Pi (a,b)>0 

-> [ne^LB < 5«] (3) 

Again this more general explanation expresses the energetic reasoning propagation and 
the bounds of activities in {i € V \ {u} \ Pi(a, b) > 0} can further be generalised in the 
same way as for a resource overload. But here the available energy units A for widening 
the bounds is rest(a, b, u) — r u ■ (newLB — 1) + 1. Hence, < A < r u indicate that the 
explanation only can further be generalised a little bit. We perform this generalisation 
as for the overload case. 



4. Experiments on Resource-constrained Project Scheduling 
Problems 

We carried out extensive experiments on Rcpsp instances comparing our solution ap- 
proach using both time-table and/or TtEf propagation. We compare the obtained 
results on the lower bounds of the makespan with the best known so far. Detailed 
results are available at http://www.cs.mu.oz.au/~pjs/rcpsp, 



We used six benchmark suites for which an overview is given in Tab. [T] where $dnst, 
#act, pi, and #res are the number of instances, number of activities, range of processing 
times, and number of resources, respectively. The first two suites are highly disjunctive, 
while the remainder are highly cumulative. 

The experiments were run on a X86-64 architecture running GNU/Linux and a In- 
tel(R) Core(TM) i7 CPU processor at 2.8GHz. The code was written in Mercury [30] 
using the G12 Constraint Programming Platform |31j . 

We model an instance as in [27] using global cumulative constraints cumulative and 
difference logic constraints (Si+pi < Sj), resp. In addition, between two activities i, j in 
disjunction, i.e., two activities which cannot concurrently run without overloading some 
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resource, the two half-reified constraints [10] b — > Si + pi < Sj and -16 — > Sj + pj < Si 
are posted where 6 is a Boolean variable. 

We run cumulative constraint propagation using different phases: 

(a) time-table consistency check in 0(n +plogp) runtime, 

(b) TtEf consistency check in 0(l 2 + n) runtime as defined in Section 3.1 



(c) time-table bounds' propagation in 0(1 ■ p + k ■ min(i?, n)) runtime, and 



(d) TtEf bounds' propagation in 0(l 2 + k-n) runtime as defined in Section 3.2 where 
k,l,n,p are the numbers of bounds' updates, unfixed activities, all activities, and 
height changes in the resource profile, resp. 

Note that in our setup phase (d) TtEf bounds' propagation does not take into account 
the bounds' changes of the phase (c) time-table bounds' propagation. For the experi- 
ments, we consider three settings of the cumulative propagator: tt executes phases (a) 
and (c), ttef(c) (a-c), and ttef (a-d). Note that phases (c) and (d) are not run if either 
phase (a) or (b) detects inconsistency. 



4.1. Upper Bound Computation 

For solving Rcpsp we use the same branch-and-bound algorithm as we used in |27j . 
but here we limit ourselves to the search heuristic HotRestart which was the most 
robust one in our previous studies [26J [27] . It executes an adapted search of [I] using 
serial scheduling generation for the first 500 choice points and, then, continues with an 
activity based search (a variant of Vsids [20 1) on the Boolean variables representing a 
lower part x < v and upper part v < x of the variable x's domain where x is either a 
start time or the makespan variable and v a value of x's initial domain. Moreover, it is 
interleaved with a geometric restart policy [35] on the number of node failures for which 
the restart base and factor are 250 failures and 2.0, respectively. The search was halted 
after 10 minutes. 

The results are given in Tab. [2] and [3} For each benchmark suite, the number of solved 
instances (#svd) is given. The column cmpr(a) shows the results on the instances solved 
by all methods, where a is the number of such instances. The left entry in that column 
is the average runtime on these instances in seconds, and the right entry is the average 
number of failures during search. The entries in column all(a) have the same meaning, 
but here all instances are considered where a is the total number of instances. For 
unsolved instances, the number of failures after 10 minutes is used. 

Table [2] shows the results on the highly disjunctive Rcpsps. As expected, the stronger 
propagation (ttef(c), ttef) reduces the search space overall in comparision to tt, but the 
average runtime is higher by a factor of about 5%-70% and 50%-100% for ttef(c) and 
ttef. Interestingly, ttef(c) and ttef solved respectively 1 and 2 more instances on J60 and 
closed the instance j 120 11 on J120 which has an optimal makespan 105. This makespan 
corresponds to the best known upper bound. However, the stronger propagation does 
not generally pay off for a Cp solver with nogood learning. 
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Table 2: UB results on highly disjunctive Rcpsps. 







J30 






J60 






#svd 


cmpr(480) 


all(480) 


#svd 


cmpr(429) 


all(480) 


tt 


480 


0.12 1074 


0.12 1074 


430 


1.82 5798 


64.25 93164 


ttef(c) 


480 


0.20 1103 


0.20 1103 


431 


2.00 4860 


64.39 80845 


ttef 


480 


0.23 991 


0.23 991 


432 


3.04 5191 


64.87 62534 






J90 






.1120 




#svd 


cmpr(400) 


all(480) 


#svd 


cmpr(280) 


all(600) 


tt 


400 


5.03 9229 


104.09 132234 


283 


9.71 15022 


322.35 398941 


ttef(c) 


400 


6.93 9512 


105.69 104297 


282 


13.47 16958 


324.73 297562 


ttef 


400 


8.10 8830 


106.66 72402 


283 


14.97 13490 


324.66 186597 






AT 












#svd 


cmpr(129) 


all(144) 








tt 


132 


8.90 19997 


66.22 87226 








ttef(c) 


130 


9.36 16466 


69.41 72056 








ttef 


129 


13.55 17239 


74.60 63554 









Table 3: UB results on highly cumulative RCPSPs. 







BL 






Pack 






#svd 


cmpr(40) 


all(40) 


#svd 


cmpr(16) 


all(55) 


tt 


40 


0.16 2568 


0.16 2568 


16 


77.65 245441 


447.69 699615 


ttef(c) 


40 


0.02 370 


0.02 370 


39 


37.22 122038 


186.79 292101 


ttef 


40 


0.02 269 


0.02 269 


39 


44.44 105751 


188.23 257747 






KSD15_ 


D 




Pack_d 




#svd 


cmpr(480) 


all(480) 


#svd 


cmpr(37) 


all(55) 


tt 


480 


0.01 26 


0.01 26 


37 


32.72 42503 


218.26 184293 


ttef(c) 


480 


0.01 26 


0.01 26 


37 


23.96 32916 


212.37 170301 


ttef 


480 


0.01 26 


0.01 26 


37 


36.93 37004 


221.11 157015 
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Table 4: LB results on AT, Pack, and Pack_d 





AT 


Pack 


Pack_d 


ttef(c) 
ttef 


5/4/3 +52 
7/2/3 +44 


0/4/12 +100 
1/4/11 +101 


0/7/11 +632 
2/5/10 +618 



Table 5: LB results on j60, j90, and J120 





J60 

+1 +2 +3 


j90 

+1 +2 +3 +4 +5 


J120 

+1 +2 +3 +4 +5 +6 +7 +8 +9 +10 


. ttef(c) 
1 mm 

ttef 


4 1 - 
7 5 - 


12 1 - - - 
25 14 3 1 - 


27 8 4 - - - 2 - - 
90 20 10 5 2 - - 2 - 


m • ttef ( c ) 
10 nuns 

ttef 


21 2 - 
13 6 3 


25 7 - - - 
35 17 6 3 1 


68 16 4 4 2 - - 1 1 
116 39 9 9 4 1 - - 1 1 



Table [3] presents the results on highly cumulative Rcpsps which clearly shows the 
benefit of TtEf propagation, especially on BL for which ttef(c) and ttef reduce the 
search space and the average runtime by a factor of 8, and Pack for which they solved 
23 instances more than tt. On Pack_d, ttef(c) is about 50% faster on average than tt 
while ttef is slightly slower on average than tt. No conclusion can be drawn on KSD15_D 
because the instances are easy for Lcc solvers. 



4.2. Lower Bound Computation 

The lower bound computation tries to solve RCPSPs in a destructive way by converging 
to the optimal makespan from below, i.e., it repeatedly proves that there exists no 
solution for current makespan considered and continues with an incremented makespan 
by 1. If a solution found then it is the optimal one. For these experiments we use the 
search heuristic HotStart as we did in [261 EI]- This heuristic is HotRestart (as 
decribed earlier) but no restart. We used the same parameters as for HotRestart. 
For the starting makespan, we choose the best known lower bounds on j60, j90, and 



J120 recorded in the PSPLib at http : //129 . 187 . 106 . 231/psplib/ and [31] at http : // 



|vilim. eu/petr7c paior2011-res ults ,txt[ On the other suites, the search starts from 
makespan 1. Due to the tighter makespan, it is expected that the TtEf propagation will 
perform better than for upper bound computation on the highly disjunctive instances. 
The search was cut off at 10 minutes as in |26| 127]. 

Table [4] shows the results on AT, Pack, and Pack_d restricted to the instances that 
none of the methods could solve using the upper bound computation, that are 12, 16, 
and 18 for AT, Pack, and Pack_d, respectively. An entry a/b/c for method x means 
that x achieved respectively a-times, 6-times and c-times a worse, the same and a better 
lower bound than tt. The entry +d is the sum of lower bounds' differences of method 
x to tt. On Pack and Pack_d, ttef(c) and ttef clearly perform better than tt. On the 
highly disjunctive instances in AT, ttef(c) and tt are almost balanced whereas tt could 
generate better lower bounds on more instances as ttef. The lower bounds' differences 
on AT are dominated by the instance stl03_4 for which ttef(c) and ttef retrieved a lower 
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bound improvement of 54 and 53 time periods with respect to tt. 

The more interesting results are presented in Tab. [5] because the best lower bounds 
are known for all the remaining open instances (48, 77, 307 in J60, J90, Jl20)Q An 
entry in a column +d shows the number of instances for that the corresponding method 
could improve the lower bound by d time periods. On these instances, we run at first 
the experiments with a runtime limit of one minute as it was done in the experiments 
for TtEf propagation in [34] but he used a Cp solver without nogood learning, tt could 
not improve any lower bound because its corresponding results are already recorded 
in the PSPLib. ttef(c) and ttef improved the lower bounds of 59 and 183 instances, 
respectively, which is about 13.7% and 42.4% of the open instances. Although, the 
experiments in were run on a slower machind^ the results confirm the importance 
of nogood learning. For the experiments with 10 minutes runtime, we excluded tt due 
to time constraints and expected inferior results to ttef(c) and ttef. With the extended 
runtime, ttef(c) and ttef could improved the lower bounds of more instance, namely 151 
and 264 instances, respectively, which is about 35.0% and 61.1%. Moreover, 3, 1, and 
1 of the remaining open instances on j60, J90, and j120, respectively, could be solved 
optimally. See App. [B]for the listing of the closed instances and the new lower bounds. 

5. Conclusion and Outlook 

We present explanations for the recently developed TtEf propagation of the global 
cumulative constraint for lazy clause generation solvers. These explanations express an 
energetic reasoning propagation which is a stronger propagation than the TtEf one. 

Our implementation of this propagator was compared to time-table propagation in 
lazy clause generation solvers on six benchmark suites. The preliminary results confirms 
the importance of energy-based reasoning on highly disjunctive RCPSPs for Cp solvers 
with nogood learning. 

Moreover, our approach with TtEf propagation was able to close one instance. It 
also improves the best known lower bounds for 264 of the remaining 432 remaining open 
instances on RCPSPs from the PSPLib. 

In the future, we want to integrate the extended edge-finding propagation into TtEf 
propagation as it was originally proposed in [33], to perform experiments on cutting and 
packing problems, and to study different variations of explanations for TtEf propaga- 
tion. Furthermore, we want to look at a more efficient implementation of the TtEf 
propagation as well as an implementation of energetic reasoning. 
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A. TTEF propagation algorithms 



Algorithm [I] shows the TtEf consistency check used. The outer loop (lines 2-16) iterates 
over all distinctive possible end times for the time windows while the inner loop (lines 
7-16) iterates over all possible start times. In line 11 (12), it checks whether a must 
fully (partially) be executed in the current time window and further ones checked in the 
same inner loop. If so it adds the required free energy units e FF of a to E. In line 13, 
it calculates the still available energy units in the time window [begin, end) taking the 
energy units from the resource profile ttEn(a, b) into account. If this results in a resource 
overload then a corresponding explanation is generated (line 15) and the algorithm fails; 
otherwise, the algorithm succeeds. 



Algorithm 1: TtEf consistency check. 

Input : X an array of activities sorted in non-decreasing order of the earliest start time. 
Input : Y an array of activities sorted in non-decreasing order of the latest completion time. 



1 end — oo; 

2 for y := n down to 1 do 

b~Y[y]; 

4 if lctf> = end then continue; 

5 end := Ictb', 
e E~0; 

7 for x :— n down to 1 do 

8 a := X[x]; 

9 if end < est a then continue; 
10 begin := est a ; 

n if lct a < end then E :— E + ef F ; 

12 if lsta F < end then E := E + r a ■ (end - Zstf F ); 

13 avail := R ■ (end — begin) — E — ttEn(a, b); 

14 if avail < then 

15 explainOverload(&egrm, end); 

16 return false; 
it return true; 



Algorithm [2] shows the lower bounds propagation algorithm. As for Alg. [T] the outer 
loop (lines 3-24) and inner loop (lines 7-24) iterate over the end and start times of the 
time windows [begin, end), but require more book keeping. In line 6, it initialises E. 
u, and enReqU where: E records the required energy units by the considered activities 
that must fully or partially be run in the time window; and u stores the activity that 
maximises min(e^, r u ■ (end — begin)) — r u ■ max(0, end — lst^ F ) and that value is saved 
in enReqU. If a must be fully or partially be executed in the time window then the 
corresponding energy units are added to E in lines 11 and 14, resp. The desired activity 
for pruning is computed in lines 13, 15, and 16, whereas the available energy units are 
calculated in line 17. In the case that there is not sufficient energy available then the 
condition of line 18 holds and the algorithm determines the first possible start time for u 
(lines 19, 20). If that is larger than the recorded earliest start time in est' u then the 
algorithm generates the explanation (line 22) and postpones the update (line 23) after 
finishing with the outer loop (line 25). 
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Algorithm 2: TtEf lower bounds propagator on the start times. 

Input : X an array of activities sorted in non-decreasing order of the earliest start time. 

Input : Y an array of activities sorted in non-decreasing order of the latest completion time. 
1 for i G V EF do est'i := esU; 



2 


end 


:= oo; k := 0; 


3 


for 


y := n down to 1 do 


4 




b:=Y[y\; 


5 




if Ictb = end then continue; 


6 




end := ic£&; i? := 0; u := — oo; enReqU := 0; 


7 




for x := n down to 1 do 


8 




a := X[x]; 


9 




if end < est a then continue; 


10 




begin := est a ; 


11 




if lct a < end then E := E + ef F ; 


12 




else 


13 




en/n := r a • max(0,end — lst EF )\ 


14 




E := E + enln; 


15 




enReqA := min(ef' F , r a ■ (end — est a )) — enln: 


16 




if enReqA > enReqU then it := a; enReqU :- 


17 




avail :— R ■ (end — begin) — E — ttEn(a, b); 


18 




if enReqU > and avail — enReqU < then 


19 




rest := E — avail — r a ■ max(0, end — lst a ); 


20 




IbU := begin + \rest/r u ~\', 


21 




if est' u < IbU then 


22 




expl := explainUpdate(begin,end,u,est' u , 


23 




Update[++k] := (u,lbU,expl); 


24 




est' u ~ IbU; 


25 


for 


z := 1 to k do updateLB(Update[z}); 



enReqA; 
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Table 6: New lower bounds on j60. 



inst LB 


inst LB 


inst LB 


inst LB 


inst LB 


inst LB 


inst LB 


9_1 85 
13_2 103 
25_2 96 
41_3 90 


9_5 81 
13_3 84 
25_4 106 
41J.0 106 


9_6 106 
13_4 98 
25_6 106 
45_1 90 


9_7 103 
13_7 82 
29_1 97 


9_8 95 
13_8 115 
29_6 144 


9_10 89 
13_9 96 
29_7 115 


13_1 105 
13_10 113 
29_8 97 



Table 7: New lower bounds on j90. 



inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


5_3 


84 


5_5 


109 


5_7 


106 


5_8 


97 


5_9 


114 


5_10 


95 


9_2 


122 


9_3 


98 


9_4 


120 


9_5 


127 


9_6 


113 


9_7 


103 


9_8 


111 


9_9 


106 


9_10 


105 


13_2 


119 


13_3 


105 


13_5 


109 


13_7 


116 


13_8 


113 


13_9 


117 


13J.0 


114 


21_7 


106 


21_8 


108 


25_1 


117 


25_2 


122 


25_3 


113 


25_4 


128 


25_5 


110 


25_6 


113 


25_8 


131 


25_9 


98 


25_10 


119 


29_1 


126 


29_2 


122 


29_4 


139 


29_6 


117 


29_7 


160 


29_8 


146 


29_9 


120 


30_9 


92 


37_2 


114 


41_1 


129 


41 _2 


154 


41_3 


149 


41_4 


142 


41 _5 


116 


41_6 


124 


41_7 


145 


41_8 


148 


41 _9 


110 


41_10 


144 


45_1 


143 


45_2 


138 


45_3 


144 


45_4 


126 


45_6 


163 


45_7 


129 


45_8 


150 


45_9 


145 


45_10 


156 


46_9 


86 







B. Closed Instances and New Lower Bounds on PSPLib 

From the open instances, we closed the instances 9_3 (100), 9_9 (99), 25_10 (108) on j60, 
5_6 (86) on j90, and 1_1 (105), 8_6 (85) on Jl20 where the number in brackets shows the 
optimal makespan. We computed new lower bounds on the remaining open instances 
from the PSPLib. Tables [6HU list these new lower bounds where the column "inst" shows 
the name of the instance and the column "LB" the corresponding new lower bound. 

C. Best Lower and Upper Bounds Retrieved 

For a later comparison, Tables |9}fTT1 show the best lower and upper bounds for AT, 
Pack, and Pack_d retrieved by one of the methods tt, ttef(c), and ttef. The column 
"inst" shows the instance name and the column "LB/UB" the corresponding lower and 
upper bound. If these bounds are equal then only one number is given. 
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Table 8: New lower bounds on j120. 



inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


inst 


LB 


6_1 


134 


6_2 


127 


6_5 


117 


6_6 


141 


6_8 


141 


6_9 


150 


6_10 


158 


7_1 


99 


7_3 


98 


7_4 


106 


7_6 


116 


7_7 


114 


7_8 


93 


7_9 


87 


7_10 


112 


8_2 


102 


8_5 


100 


8_9 


90 


8_10 


92 


9_4 


85 


11 1 


157 


11_2 


147 


11_3 


189 


11_4 


178 


11_5 


194 


11_6 


192 


11_7 


149 


11_8 


153 


11_10 


164 


12_1 


126 


12_2 


112 


12_4 


122 


12_5 


155 


12_6 


116 


13_1 


124 


13_3 


116 


13_4 


109 


13_6 


96 


13_9 


83 


14_2 


91 


14_5 


94 


14_7 


90 


16 1 


181 


16_3 


221 


16_4 


191 


16_6 


195 


16_8 


183 


17_5 


124 


17_6 


134 


18_8 


102 


18_9 


89 


18_10 


97 


26_1 


155 


26_2 


159 


26_3 


158 


26_4 


161 




1 on 

ioy 


Of; R 

ZO_0 


1 71 
1(1 


Of; 7 


1/17 


Of? o 


1 fin 

loo 




1 P. 1 


of; i r\ 
ZO-W 


1 *7Q 

1 to 


07 1 


1 H7 

1U ( 


27_2 


110 


27_3 


142 


27_4 


105 


27_5 


106 


27_6 


133 


27.7 


119 


27_8 


136 


27_9 


121 


27_10 


111 


28_1 


106 


31_1 


181 


31_2 


176 


31_3 


160 


31_4 


195 


31_5 


187 


31_6 


182 


31_7 


191 


31_8 


176 


31_9 


176 


31_10 


202 


32_1 


144 


32_2 


123 


32_5 


133 


32_6 


122 


32_8 


132 


33_1 


105 


33_2 


107 


33_3 


102 


33_4 


107 


33_8 


107 


33_9 


109 


34_1 


76 


34_2 


103 


34_3 


99 


34_5 


102 


36_1 


201 


36_3 


218 


36_5 


213 


36_7 


196 


36_9 


203 


37_2 


141 


37_5 


195 


37_8 


169 


37_9 


138 


38_1 


105 


38_2 


119 


38_4 


138 


38_6 


119 


38_7 


103 


38_10 


137 


39_2 


105 


40_1 


80 


42_1 


107 


46_1 


172 


46_2 


187 


46_3 


163 


46_5 


136 


46_7 


158 


46_9 


157 


46_10 


175 


47_1 


130 


47_3 


119 


47_4 


120 


47_5 


126 


47_6 


128 


47_7 


114 


47_8 


124 


47_10 


128 


48_4 


123 


51 1 


186 


51_2 


200 


51_3 


193 


51_4 


197 


51_6 


193 


51_7 


185 


51_8 


186 


51_9 


190 


51.10 


201 


52_1 


161 


52_2 


169 


52_3 


126 


52_4 


157 


52_5 


158 


52_6 


183 


52_7 


142 


52_8 


148 


52_9 


142 


52_10 


131 


53_1 


138 


53_2 


109 


53_4 


138 


53_5 


109 


53_6 


101 


53_8 


135 


53_10 


124 


54_1 


102 


54_5 


107 


54_6 


104 


54_8 


100 


54_9 


105 


57_1 


173 


57_2 


151 


57_3 


176 


57_5 


170 


57_6 


176 


57_7 


156 


57_9 


157 


58_2 


122 


58_3 


117 


58_4 


138 


58_5 


116 


58_6 


135 


58_7 


143 


58_8 


126 


58_9 


126 


59_5 


104 


59_6 


112 


59_8 


107 


59_9 


117 


59_10 


128 


60_3 


88 


60_7 


91 
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Table 9: Lower and upper bounds for AT. 



inst 


LB/UB 


inst 


LB/UB 


inst 


LB/UB 


inst 


LB/UB 


inst 


LB/UB 


27.1 


41 


27_2 


53 


27_3 


68 


27_4 


112/114 


27.5 


56 


27_6 


73 


27_7 


54 


27_8 


95 


27_9 


38 


27.10 


45 


27_11 


57 


27.12 


73 


27.13 


38 


27.14 


55 


27.15 


46 


27.16 


75 


27.17 


55 


27.18 


55 


27.19 


79 


27.20 


152 


27_21 


92 


27_22 


86 


27_23 


82 


27_24 


106 


27.25 


51 


27_26 


53 


27_27 


58 


27_28 


95 


27_29 


51 


27.30 


76 


27_31 


75 


27_32 


82 


27_33 


66 


27_34 


61 


27.35 


115 


27.36 


146 


27.37 


78 


27.38 


100 


27.39 


119 


27.40 


130 


27.41 


60 


27.42 


53 


27.43 


75 


27.44 


88 


27.45 


49 


27_46 


65 


27_47 


75 


27_48 


80 


51_1 


98 


51.2 


96 


51_3 


133 


51_4 


161/219 


51_5 


97 


51_6 


126 


51.7 


120 


51_8 


194 


51_9 


74 


51_10 


73 


51.11 


99 


51.12 


116/137 


51_13 


84 


51_14 


86 


51_15 


86 


51.16 


132 


51.17 


84 


51_18 


99 


51_19 


170 


51.20 


274 


51.21 


145 


51.22 


168 


51 23 


183 


51.24 


228 


51.25 


95 


51.26 


89 


51.27 


113 


51_28 


164 


51_29 


98 


51_30 


105 


51.31 


130 


51.32 


139 


51_33 


116 


51_34 


115 


51_35 


173 


51.36 


300 


51.37 


162 


51_38 


177 


51_39 


189 


51_40 


218 


51.41 


102 


51.42 


108 


51.43 


121 


51.44 


174 


51.45 


122 


51.46 


125 


51.47 


151 


51.48 


167 


103 1 


158 


103.2 


182 


103.3 


216/259 


103.4 


280/445 


103_5 


191 


103_6 


207/209 


103_7 


234/293 


103.8 


207/294 


103.9 


139 


103_10 


119 


103.11 


160/169 


103_12 213/302 


103.13 


127 


103.14 


152 


103_15 157/168 


103_16 


167/179 


103_17 


209 


103.18 


232 


103.19 


301 


103.20 


475 


103.21 


276 


103.22 


295 


103.23 


368 


103.24 


449 


103.25 


177 


103.26 


183 


103.27 


199 


103.28 


295 


103.29 


225 


103_30 


231 


103_31 


227 


103_32 


281 


103.33 


220 


103.34 


264 


103_35 


341 


103_36 


575 


103_37 


327 


103.38 


376 


103.39 


389 


103_40 


451 


103_41 


191 


103_42 


187 


103.43 


260 


103.44 


375 


103_45 


216 


103_46 


251 


103_47 


262 


103.48 


300 







Table 10: Lower and upper bounds for PACK. 



inst 


LB/UB 


inst LB/UB 


inst LB/UB 


inst 


LB/UB 


inst 


LB/UB 


inst 


LB/UB 


001 


23 


002 


32 


003 


29 


004 


43/44 


005 


42 


006 


47 


007 


41 


008 


44 


009 


57/72 


010 


38 


011 


44 


012 


45 


013 


36 


014 


45 


015 


43 


016 


63 


017 


62 


018 


60 


019 


59 


020 


62 


021 


51 


022 


59 


023 


51 


024 


56 


025 


69/70 


026 


54 


027 


55 


028 


64 


029 


43 


030 


20 


031 


70 


032 


80 


033 


78 


034 


73 


035 


73/77 


036 


100/106 


037 


116/138 


038 


86 


039 


99/111 


040 


87/91 


041 


27 


042 


29 


043 


105 


044 


103 


045 


86/87 


046 


110/128 


047 


103/107 


048 


76/77 


049 


29 


050 


94/109 


051 


29 


052 


85 


053 


97/113 


054 


92/100 


055 


91/97 























21 



Table 11: Lower and upper bounds for Pack_d. 



inst 


LB/UB 


inst LB/UB 


inst 


LB/UB 


inst 


LB/UB 


inst 


LB/UB 


001 


612 


002 745/747 


003 


624/625 


004 


1381 


005 


983 


006 


1119 


007 1082 


008 


1274 


009 


1593/1951 


010 


1216 


Oil 


940 


012 1234/1241 


013 


829 


014 


1565 


015 


1198 


016 


1783/1813 


017 1641/1651 


018 


1462/1480 


019 


1526/1542 


020 


1661 


021 


1606 


022 1787 


023 


1092 


024 


1625 


025 


2061/2147 


026 


926 


027 1789/1793 


028 


1897/1962 


029 


1233 


030 


597 


031 


1949 


032 2943 


033 


3390 


034 


2371 


035 


2305 


036 


2175/2191 


037 3325/3614 


038 


2180 


039 


2730/2734 


040 


3024 


041 


679 


042 838 


043 


2439 


044 


3050 


045 


2712 


046 


3243/3277 


047 2740/2745 


048 


2446 


049 


675 


050 


2687/2716 


051 


838 


052 2253 


053 


2521 


054 


2750 


055 


2628 
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